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Abstract 

The HNC Ornstein-Zernike integral equations are used to determine the properties of simple 
models of colloidal solutions where the colloids and ions are immersed in a solvent considered as 
a dielectric continuum and have a size ratio equal to 80 and a charge ratio varying between 1 
and 4000. At an infinite dilution of colloids, the effective interactions between colloids and ions 
are determined for ionic concentrations ranging from 0.001 to 0.1 mol/1 and compared to those 
derived from the Poisson-Boltzmann theory. At finite concentrations, we discuss on the basis of 
the HNC results the possibility of an unambiguous definition of the effective interactions between 
the colloidal molecules. 
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I. INTRODUCTION 



The large difference of sizes and electric charges between the colloids and molecules or 
atoms present in the colloidal suspensions is one of the main difficulties encountered in 
the theoretical study of these multicomponent solutions. A customary simplification is 
to approximate the solvent contributions to the solution properties by those of a dielectric 
continuum where are immersed the colloids and ions. From this simplification, it results that 
the colloidal suspensions are mixtures of colloids and ions where the coulombic interactions 
between charges are divided by the dielectric constant of the continuum. All excluded volume 
and polarization effects due to the solvent are neglected or assumed to increase notably the 
ion sizes by considering that the ions are surrounded with a shell of bound molecules of 
solvent, for instance a hydration shell. A new simplification is generally done since the 
colloids are considered as structureless, rigid, spherical and charged particles. With these 
approximations, the estimate of properties of colloidal suspensions is reduced to that of 
mixtures of spherical ions and particles where the size and charge of ions and particles have 
highly different values |l[ O, . 

The screening effects which are present in the coulombic systems, lead to further simplify 
the description of the colloidal suspension by supposing that the interaction between colloidal 
particles due to their electric charges can be represented by an effective potential resulting 
from the screening of their own charges by the ions 

mm 

. The standard methods of 

liquid theory allow to avoid this latter approximation and to calculate the structural order 
and thermodynamic properties of the colloidal suspensions as those of mixtures of large 
charged particles and small ions. The main aim of this work is to demonstrate that the 
HNC Ornstein-Zernike (OZ) integral equations j3| allow to perform such a calculation for 
a size ratio between spherical colloidal particles and spherical ions equal to 80 and charge 
ratios ranging from 1 to 4000. Although these charge ratios are currently found in real 
suspensions, the size ratio is small compared to that characteristic of real colloids. However, 
a size ratio equal to 80 preserves the fact that the external surface of colloids is sufficiently 
large for allowing that the screening of the colloid charge by the ions is not hindered by 
excluded volume effects, because the external colloid surface is able to accommodate about 
25000 small ions. Such a size ratio seems sufficient to perform a valuable comparison with the 
Poisson-Boltzmann (PB) theoretical approaches of the estimate of the effective interaction 
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between ions and colloids in the infinite dilution limit, in particular for values of the charge 
ratio larger than 1000. 

In addition to the coulomb interactions, short range interactions of van der Waals type 
exist between the ions and colloids. These interactions will be also discarded in order 
to determine unambiguously the effective interactions induced solely by the interactions 
between the charges. The specific system studied is a neutral mixture of large charged hard 
spheres, the colloids, soluted in a fluid of one or two species of small charged hard spheres, 
the ions. The colloids are positively charged. A link between this model system and the 
real suspensions is made by attributing to the ions a hard core diameter of cxj = 5 A, a 
typical value for hydrated ions, giving for the colloid diameter a value of cTc = 400 A; this 
equivalence allows to express the ionic densities of the model in mol/1 (M). 

The numerical procedure used to solve the HNC integral equations is described in Section 
II, in the same section the results obtained at an infinite dilution of colloids are given and 
compared to those of the PB theory. In Section III we present the results obtained at finite 
concentration of the colloids. In the final Section, we summarize the results and discuss 
possible extensions of this work. 



II. INFINITE DILUTION 



To solve numerically the HNC OZ equations for symmetric and asymmetric mixtures of 
charged hard spheres, a sufficiently small integration step is needed to describe accurately the 
very steejp peak existing at contact in the correlation functions between the negative ions and 
co,loid.fi6QaflyQ.Fu.he..o.e..Keuppe. bound of ..e„Hn.e..a.o. 
must be equal to several colloid diameters in order to study systems with a colloid finite 
concentration. Finally, at low ionic concentrations, the Debye length characterizing the 
screening effects can have a value of the same order than the colloid diameter which requires 
also the use of an upper bound of integration large enough to reach the asymptotic behavior 
of the correlation functions. With these constraints and the fact that the numerical solution 
of the HNC equations is obtained using a Fast Fourier transform algorithm, the numerical 
integrations are performed with 2^^ = 262144 Ar or Ak steps in the r and k spaces. The 
colloid radius Rc, chosen as the unit of length, is equal to 16000 Ar and the ion radius to 
200 Ar. In r space, the integration upper bound is equal to 16.38 Rc- 
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FIG. 1: Ion-ion correlation functions gn{r), guir) (solid lines) and g22{r) (filled squares) with 
zi = —1, Z2 = 1, pv = 0.01 M and po = 0. Insert : fit of gii{r) and 512 (r) (solid lines) by 
Aaf)e^p{—Kpr)/r (filled circles) {KpRc = 6.62). 

The coupling strength between charges is characterized by the Bjerrum length Ib — 
e^/ekBT where e is the dielectric constant of the solvent, e the electron charge, ks the 
Boltzmann constant and T the temperature. For an aqueous suspension Ib — 7.198 A at 
T = 293 K. The colloid and ion densities pO) Pi and p2 must be chosen in order to satisfy 
the electro-neutrality relation zq po + zipi + Z2P2 = where zq, z\ and Z2 are the colloid 
and ion charges. 

In the solution of the HNC equations, the long range of the coulomb potential is taken 
into account by writing the direct correlation functions Capir) (a, /3 = 0, 1 and 2) : 

ui{air) 



Ca/sir) = aa/3{ai, r) - Ib z^ zp 



(1) 



where ui{air) — erf (a^r) and a; is chosen equal to 0.9 Re- The functions ao,p{ai,r) have a 
short range, because the asymptotic behavior of 0^/3 (r) at large r is equal to —Ib ZaZ/j/r. 
The Fourier transform f{k) of a function /(r) being defined by 

(2) 



kr 

where p = po + Pi + P2 (for an infinite dilution of colloids po = 0), the Fourier transform of 
ui{air) is ui{k) — exp{—k'^/{4af))/k'^ and that of CQ,^(r) 

Cafs{k) = aa/3{k) - ciui{k) z^zp. (3) 




FIG. 2: Ion-colloid correlation functions goi{r) (solid line) and go2{r) (dashed line) with zi = — 1, 
Z2 = 1, Pv = 0.001 M and po = 0. Insert : fit of goiir) (solid line) and go2{r) (dashed line) by 
Aaf3ex.p{—Kpr)/r (filled circles) (KpRc = 2.08). 



where q = An pis- 

In /c-space, the OZ equations at a colloidal infinite dilution write for the ion correlation 
functions 

Cii{k) + X2{cl2{k) - Cu{k) C22{k)) 



hu{k) 
h22{k) 
Jiuik) 



Do{k) 

C22{k) + Xi(Ci2(/c) - Cii{k) C22{k)) 



Do{k) 



cuik) 



Do{k) 

and the ion-colloid correlation functions 

Coi{k) (1 -X2C22ik)) +X2 Ci2{k)co2{k) 



(4) 



hoi{k) 
ho2{k) = 



Do{k) 

Co2{k) (1 - xi cu{k)) + xi Ci2{k) Coi{k)) 



Do{k) 

where Xi — pi/ p, X2 — P2/ P and 



(5) 



Do{k) = (1 - xi cii(/c)) (1 - X2 C22{k)) - xi X2 Ci2(A:) 



(6) 
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FIG. 3: At pq = 0, PB and HNC estimates of Z^ff versus zq for ionic charges zi = —1, 2:2 = 1 at 
p„ = 0.001 M : PB dotted line, HNC solid line and filled diamonds, = 0.01 M : PB dashed line, 
HNC solid line and filled circles and = 0.1 M : PB dash-dotted line, HNC solid line and filled 
squares. 

supplemented by the HNC closure in r-space : 

hapir) = Qafsir) - 1 

= exp{-lBZa zp/r + Kp{r) - Cap{r)) - 1 
= exp(-/BZ„ zpil - erf(azr))/r) 

X exp(/ia/3(r) - aa/3(ai,r)) - 1 . (7) 

The convergence of the numerical solutions is obtained by using a Picard algorithm and 
considered to be achieved when the two successive iterations n and + 1 of the functions 
lapir) = hafsir) — aa/3{ai,r) satisfy the relations: 

l7i";^^(r) - 7l"j(r)| < 10"^ (8) 

In the literature, detailed comparisons of the two-body correlation functions and ther- 
modynamic properties computed from the HNC theor y a nd Monte-Carlo simulations have 



been made for ionic solutions 
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17 



19 



They have defined the do- 



main of validity of the HNC approximation for ionic concentrations ~ 1 M and charge and 
size ratios of the order of 1-20. These works have shown that the HNC closure gives an 
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FIG. 4: At pq = 0, PB and HNC estimates of Z^jf versus zq at = 0.001 M, ionic charges zi = 
— 1, Z2 = 2, : PB solid line and white circles, HNC solid line and filled circles, and zi = —2, Z2 = 1 
: PB solid line and white squares, HNC solid line and filled squares. 

accurate description of the correlations at large r in symmetric or asymmetric systems of 
charged hard spheres, with equal or different diameters, immersed in a dielectric continuum. 
This satisfactory result constitutes a favorable basis to the study of the infinite dilution of 
colloids in our model. At po = 0, the HNC equations ^ for the correlation functions of 
the small ions are independent of those giving the ion-colloid correlation functions goi{f) 
and go2{r) (cf. Eqs. (0))- Hence, it is possible to determine goi{r) and go2{r) at po = 
for any charge zq when the equations Q have been solved. Such a computation has been 
realized for zq < 4000 and symmetric {zi = —1, Z2 = 1) or asymmetric {zi = —2, Z2 = I) 
and {zi = —1, Z2 = 2) ionic solutions where the densities of the ionic species with the 
highest valence, are equal to 0.1, 0.01 and 0.001 M (i.e. p^af = 0.007528, 0.0007528 and 
0.00007528), respectively. 

From the correlation functions, there are several possible definitions of the effective in- 
teractions between colloids or more generally between specified components of a multi- 
component mixture One of those is the potential of mean force v^^Jr) given by 
\og{gai3{r)} = —(3v'^p{r). Theoretical approaches, such as the DLVO theory |l|,y|, provide 
an explicit form for the effective interaction between colloids : 



where is the inverse Debye length defined by k|, = AnlBizfpi + Z2P2), and the parameter 




(9) 
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FIG. 5: At pq = 0, = 0.001 M and ionic charges zi = —1, Z2 = 2, variation of the functions 
goiir) near r R^: zq = 1200 sohd line and filled squares, zq = 1800 solid line and filled diamonds, 
zq = 2400 solid line and filled circles. Insert : for l.bRc < r < 3Rc, goi{f) (solid line and filled 
squares zq = 1200, diamonds zq = 1800 and circles zq = 2400) and go2{r) (dashed line and filled 
squares zq = 1200, diamonds zq = 1800 and circles zq = 2400). 

Zeff is conventionally called the effective charge of the colloids. 

In the limit po = 0, from the linearized PB theory, a similar expression is obtained for 
the effective interaction between colloid and ions 
exp(Kz) Rc) exp(-KDr) 



eff, 
Oo I 



ZeffZa Ib 



(10) 



^1 + kdRc) r 

where a = 1 or 2. For symmetric and asymmetric ionic solutions, without colloids, the PB 
theory also predicts that, at large r, the effective interactions between the ions are Yukawa 
potentials similar to fg{/(r). This result, for instance, has been verified by computing f ™^(r) 
from the correlation functions gap{f) determined by MC simulations for systems of charged 
hard spheres with equal diameters and densities of ~ IM 22] . Similarly, an analytical 
expression of v'^p{r), with the same functional dependence at large r as fo{/, will be found 
from the correlation functions computed from the HNC OZ equations Q and (0), when 
the singularity nearest to the real axis of the function l/Z)o(^) is a simple pole located at 



k 
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Numerical methods have been derived for the determination of the zeros of -Do(fc) for 



complex values of k 
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They need the calculation of integrals having the 
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FIG. 6: At po = 0, = 0.001 M and ionic charges zi = —2, Z2 = 1, variation of the functions 
goi{r) near r Rc zq = 1200 solid line and filled circles, zo = 1800 solid line and filled squares, 
zo = 2400 solid line and filled diamonds. Insert : for 1.5i?c < r < 3Rc, goi{r) (solid line and filled 
circles zq = 1200, squares zq = 1800 and diamonds zq = 2400) and 502 (?^) (dashed line and filled 
circles zq = 1200, squares zq = 1800 and diamonds zq = 2400). 

typical expression 



In practice, such integrals are not easily estimated when the upper limit of the numerical 
integration on r is ~ 500 times the radius of the ions. 

This calculation requires at least an accuracy of the order of 10^*^ — 10^^'^ on the c^Jjsir) 
functions at large r which is not achieved with the present scheme used to solve the OZ 
equations. Therefore, we analyse the asymptotic behaviors of the 5'a,/3(^) at infinite dilution 
by supposing that the pole of l/Do{k), closest to the real axis, is effectively on the imaginary 
axis. With this hypothesis, Pv'^p{r) should be equal to Ao,pex.^{—Kpr)/r with Kp identical 
for all the functions v'^p{r). The determination of and Kp is obtained from a fit of 
log{(7Q,/3(r)} by exp(— ft;pr)/r, for r sufficiently large. However, the domains of large r 
values where the functions | \og{gaf3(r)}\ are smaller than 10~^ must be excluded due to the 
numerical uncertainties on the solutions of the OZ equations. The domains of r values, used 




where (fj^isir) — Ca,/3{r) + IsZaZb/T and ki and ko are the real and imaginary parts of k. 
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in the fits, correspond to the distances for which 10~^ < | \og{ga i3{r)} \ < 10~^ where the 
numerical uncertainties are two orders of magnitude smaller than the values of the fitted 
functions. 

Examples of curve fittings of \og{gai3{r)} are given in Figs. 1 and 2 ; these figures show 
unambiguously that these functions are accurately described by the presupposed analytical 
form Ac 13 exp{— Kpr) / r . For a given ionic density, the values of Hp obtained from the curve 
fitting of the functions log{5'oi('^)} and \og{go2{r)} are independent of zq when the values of 
Zq vary from 1 to 4000 and identical to those determined by the fit of log{gii{r)}, \og{gi2{r)} 
and \og{g22{r)}. At = 0.1, 0.01 and 0.001 and ionic charges {zi = —1, Z2 = 1), {zi = 
—2, Z2 = 1) and {zi = —1,-22 = 2), the values for Hp differ by a few percent from Kd 
computed for the same ionic concentrations and charges. At these weak packing fractions 
of the ions, the differences between Kp and kd result from contributions associated with the 
finite ionic density, Kp being expected to be equal to only in the limit of very low densities 
where the functions c^'^^(r) go to zero. 

In the HNC theory, the description of goi{r) and go2{f') at large r by a Yukawa function 
is in agreement with the expressions of the PB theory for the effective interactions between 
ions and colloids at infinite dilution. It is thus possible to compare the values of the pa- 
rameter Zeff recently computed from PB theory with those deduced from the values of the 
parameters Aqi and Aq2 by writing Aq^ in a form identical to that derived from PB theory 

Obviously, this expression does not result from a boundary condition similar to that 
satisfied by the solution of the PB equation on the colloid surface and can be considered as 
a definition of Z^^^^ . The possibility to write Aq^ in such a form is partly confirmed by the 
fact that Aqi/Aq2 = z\j Z2 within the limits of the accuracy of the fit. 

Recently PB expressions of Z^^^ have been derived. The most simple expression jsij, 
valid in the hmit Zq and Rc/cn — cxd, is 
4i? 

Z:^^ = -^{l + Kj,R^). (13) 
Lb 

An expression, a priori valid for all values of Zq and Rc when, at infinite dilution, the 
colloids are immersed in low density ionic solutions of monovalent ions, has been derived 
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FIG. 7: At pq = 0, = 0.001 M and ionic charges zi = —2, Z2 = 1, variation of the functions 
go2{r) near r Rc : zq = 1200 dash-dotted line, zq = 1800 dashed line, zq = 2400 solid line. 



32j | by using the asymptotically exact solution of the PB equation 



33| 



'eff 



■t(x) + 



(14) 



where x = zqIb/{Rc {'^^kdRc + 2)) and t{x) = (a/ (1 + x'^) — l)/x. This derivation has been 
recently extended to the cases of ion mixtures with asymmetric charges zi = —2, Z2 = 1 
and zi = —1, Z2 = 2 [s^. It gives Z^jj- = Zof{x) where f{x) is an intricate function of 
X = zqIb/ i^D R^i depending on the charge asymmetry of the ions. The main results of these 
theoretical estimates are that, for zq = 1 — 100, Zg// ^ Zq and, for zq > 1000, Z^ff goes to 
a constant value of the order of Z^j-jr. 

From the fits of the functions \og{gai3{r)}, similar in accuracy to those plotted in Figs. 1 
and 2, we have computed the values of Hp and Z^jj'~^ resulting from the HNC OZ theory (cf. 
Eq. (HH)). At = 0.001 M, i^pRc and i^dRc are respectively equal to 2.089 and 2.087, at 
pv = 0.01 M equal to 6.62 and 6.60, and at = 0.1 they differ by 5% and are equal to 21.8 
and 20.8. Fig. 2 shows that for r > 2Rc, goi{r) and (7o2('") are symmetric with respect to 1, 
and Fig. 3 that the PB and HNC estimates Z^jrjr and Z^^^^ are in good agreement at the 
ionic concentrations of 0.001 and 0.01 M when zq < 2000, in particular the two estimates 
reach a maximal value near zq ^ 2000. For larger values of Zq the PB and HNC estimates 
ZKf and Z^/f'-^ differ systematically : ZK. stays constant whereas Z^/f'-^ decreases. 



At pt, = 0.1 M, in the considered range of zq, the PB values of Z^ff do not reach their 
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FIG. 8: At pq = 0, PB and HNC estimates of Z(.ff versus zq for ionic concentrations = 0.01 M, 
ionic charges = — 1, 2:2 = 2, : solid line and white circles PB, solid line and filled circles HNC, 
and zi = — 2, 2:2 = 1 : solid line and white squares PB, solid line and filled squares HNC. 

asymptotic value. The HNC estimates of Z^ff would have values exceeding largely zq if 
'^e//^'" was calculated from Eq. (fT^ by using and not Kp. This result seems to be 
an indication that the finite concentration effects of co-ions and counter-ions cannot be 
neglected when the ionic concentrations are of the order of or larger than 0.1 M. 

For the charge asymmetric ionic mixtures of = 0.001 M, Zi = —2, Z2 = 1 and Zi = 

— 1, Z2 = 2, Fig. 4 presents a comparison between Z^^j and Z^^^'~^ . For the case Zi = 

— 1, Z2 = 2, the qualitative agreement is excellent and the quantitative difference is not 
larger than 10%. The two estimates of Z^jf are almost constant when zq > 1000. Fig. 5 
shows that, in agreement with the fact that Z^^^^ is quasi constant for these values of zq, 
goi{r) and go2{i^) vary only near r c:^ Rc and stay almost unchanged when r > 2 Rc. In 
the case zi = —2, Z2 = 1, Z^jj- and Zfj^^ differ for zq > 500, Z^j^j is constant but Z^j^^'-^ 
decreases and becomes negative near zq = 1800. Fig. 6 shows the variation of goi{r) and 
902(1^) which, taking into account the ratio Zi/z2 and the sign inversion of Z^j^^ , retain the 
correct symmetry with respect to 1 when r > 2 

For r < 2 /2c, the main feature is the appearance of a peak in go2{r) (cf. Fig. 7) near 
r ~ i?c + l.ScTj characterizing an overscreening of the colloid charge for large zq by the 
negative ions which allows the positive ions to be located close to the colloid surface. This 
overscreening seems to be at the origin the sign inversion of the colloid-ion potential of mean 



12 



o 




_i I I I I I I i_ 



1.0 1.1 1.2 



1.00 1.02 1.04 1.06 1.08 1.10 
r/R 



FIG. 9: At po = 0, pt, = 0.1 M zi = —2, Z2 = 1 and zq = 3800 : ion-colloid correlation functions 
goi{r) and 5o2(?")- 

force, which is not obtained in the PB theory. 

A similar behavior is found for = 0.01 as shown by the comparison between Z^jj and 
Zg^y^ presented in Fig. 8. At p„ = 0.1, zi = -2, Z2 = I {KpR^ > 25), Fig. 9 shows that 
5'02(^) exhibits a peak at r ~ + l-5crj for zq = 3800 which indicates an overscreening 
of the colloid charge. However the large value of UpRc implies a very steep decrease of the 
correlation functions and the range of r values where a significant fit of goi{r) and go2{f') is 
possible, is very small O.lRc) and does not allow for an estimate of Z^j^^'^ . 

III. FINITE CONCENTRATION 



The binary mixture composed of colloids and counterions with zi = —1, without added 
salt, has been well studied in the literature 0,1^ and we briefly present and discuss the HNC 
results obtained for the case of a coUoid-counterion size ratio equal to 80 at a colloid packing 
fraction rj^ = TTpcCr^/Q = 0.1. For this mixture, a numerical solution of the HNC OZ equations 
is only found when zq < 130. This limitation on the Zq values is similar to that obtained in 
Refs. 0, At Zq = 130, using the previously adopted values of ion and colloid radii, 

the ionic concentration is ^ 0.0064 M and the ionic packing fraction ^ 0.0000253. Even 
at low values of zq the asymptotic behavior of the three correlation functions corresponds 



13 



TABLE I: Colloid packing fraction : r/c, density of the ion species with the highest valence : 
ion charges : zi and Z2, maximum colloid charge : z™^^, k/j and : inverse screening Debye 
length without and with the colloid charge contribution. The last line corresponds to a binary 
ion-colloid mixture. 
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to damped oscillations (cf. Fig. 10) indicating that, in the OZ equations (cf. Eq. (jlj)), 
the asymptotic behavior is determined by poles of l/Do{k) symmetric with respect to the 
imaginary axis. A fit of the pair correlation functions at zq = 130 for 8Rc < r < 25Rc 
(cf. above and below for a detailed discussion of the fitting procedure) shows that these 
functions oscillate with the same period equal to 1.95i?c (cf. Fig 11) and have the same 
damping constant UpRc = 0.24 about 6 times smaller than the inverse screening Debye 
length kd Rc = 1-67 (cf. Table 1). 

In Fig. 10, it is seen from goi{r) that the screening of the colloid charges occurs over a 
distance which varies between 1.5 and 0.6 Rc when zq increases from 1 to 130. For these 
values of zq, the range of the screening effects stays sufficiently large so that the coulombic 
repulsion shifts the main peak of goo{r) from 2Rc to 3.8-Rc- This displacement induces a 
reinforcement of the excluded volume effects between the colloids. For — 50 — 130, 5'oo('") 
has the typical form of the pair correlation function of a dense monoatomic fluid. It seems, 
although difficult to establish numerically, that this increase of the excluded volume effects 
is the main cause of the instability of the HNC OZ equations at Zq > 130. 

The ternary solutions of colloids, positive and negative ions have been considered for 
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FIG. 10: At 77c = 0.1, p2 = and zi = —1 colloid-colloid goQ{r) (solid line), ion-colloid gm{r) 
(dashed line) and ion-ion gii{r) (dash-dotted line) correlation functions, from bottom to top zq = 
10, zq = 50 and zq = 130. 

different colloid packing fractions, ionic concentrations and colloid charges summarized in 
Table 1. In these mixtures for given densities of the colloids and the positive ions, the 
density of the negative ions is fixed by the electroneutrality condition. In Table 1, z^"'^ 
denotes the maximum value of Zq at which the HNC equations can be solved by using the 
previously described numerical iterative scheme. In the mixtures, the asymptotic behaviors 
of the correlation functions Qaiiij) can be analysed similarly to the infinite dilution case, 
since a typical OZ equation for a correlation function has the form 

(coi(A;) + X2 C02{k) Ci2(fc) - X2 C22{k) Coi{k)) 



D{k) 



(15) 



15 



1.00005 



^ 1.00000 

00 



0.99995 



TT 



— 1 — 
f\ 
I \ 
I \ 
I \ 



1.005 



^ 1.000 



0.995 



.15 



20 



25 



30 



/ 

I / 
1 / 
\ I 

L-L. 




1.005 



^ 1.000 

00 



0.995 




FIG. 11: At rjc = 0.1, /O2 = 0, zi = —1 and zq = 130, : colloid-colloid goo{r) (solid line), ion-colloid 
goiif) (dashed line) and ion-ion gii{r) (dash-dotted line) correlation functions. From bottom to 
top : goo{r), goi{r) and (711 (r), and, after an appropriated shift of goi{r) and gu{r) on r equal 
to (Iqi and cl^2 I^l- El)> for 4 < r/Rc < 20 and 15 < r/Rc < 30 close views showing the 
identity of the oscillation period of the three correlation functions (the maxima and minima of the 
oscillations of 500 (^) and goi{r) can be truncated). 
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FIG. 12: At r]c = 0.001, zi = -1, Z2 = I, Pv = 0.001 M and zq = 300, bottom : correlation 
functions gn{r) (solid line), gi2{r) (dashed line) and g22{r) (dash-dotted line), insert : close view 
for r > ARc, and top : goo{f) (solid line), goi{r) (dashed line) and go2{r) (dash-dotted line), insert 
: close view for r > 6Rc. 

where the expression of D{k), identical for all ha/3{k), is 
D{k) = {1 - X2C22{k) - xicii{k) - xoCoo{k) 

+ Xi X2 Cu{k) C22{k) — X1X2 cuik)"^ 

- XoXiCoi{k)'^ - XoX2Co2{k)^ 

+ XoXiCii{k)coo{k) - xo Xi X2 Cu{k) Cqq {k) C22 {k) 
+ xoXiX2Cii{k) co2(A;)^ 

- 2 Xo Xi X2 Coi{k) Ci2{k) Co2{k) 

+ XoX2Coo{k)c22{k) + XoXiX2Coi{k)'^ C22{k) 
+ XoXiX2Ci2{ky Coo{k)) 



(16) 



where Xq — po/p- 
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FIG. 13: At r]c = 0.001, zi = -2, Z2 = I, py = 0.001 M and zq = 259, bottom : correlation 
functions gii{r) (solid line), guir) (dashed line) and g22{r) (dash-dotted line), insert : close view 
for r > 4i?c, and top : goo{f) (solid line), goi{r) (dashed line) and go2{r) (dash-dotted line), insert 
: close view for r > SRc (the scale of y-axis precludes to identify the maxima and minima for 
r > GRc). 

Clearly, the zeros oi D{k), if they exist at complex values of k near the real axis, determine 
the decrease of all the functions Qaisir). Since these functions are real, the most simple 
locations of these zeros correspond to values of k on the imaginary axis, k = iKp, giving 
at large r an exponential decrease of gapif) — 1 + Ao,pex^{—Kpr)/r, or to values of k 
symmetric with respect to the imaginary axis, k — kQ-\- iKp and k — —ko + iKp, giving an 
oscillatory exponential damped decrease of gap{f) — l + ^a/3exp(— Kpr) cos{kQr-\-dap)/r. In 
the first case, at large r all correlation functions are exponentially decreasing, in the second 
case all the functions gapif^) oscillate with an identical period and damping. Contrary, to 
the infinite dilution case where, for the considered ionic solutions, an exponential decrease 
was always obtained, both exponential and oscillatory damped decreases are found for finite 
colloid concentrations. 
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For instance, in the case of a symmetric ionic solution {zi — —1, Z2 — 1), at pv — O.OOIM 
and ric — 0.001, the OZ equations can be solved in the domain 1 < 2^0 < 300. For zq < 100, 
goo{r) increases monotonically from 0, at short distance, towards 1 at large r. Above this 
value of Zq (cf. Fig. 12), goo{r) presents, at r ~ 6 — 7Rc, a weak broad peak, larger 
than 1, then, at r ~ lORc, a minimum smaller than 1, and finally, at large r, increases 
monotically towards 1. A similar behavior is obtained for gii{r) ; the peak and minimum, 
which exist also for zq > 100, are located at r ~ 1 — 2Rc and at r ~ 5 — 6Rc, respectively. 
At larger r, gii{r) increases monotically towards 1. The function g2 2{T^), different from 
5'ii(r) due to the finite density of the colloids, presents a peak and minimum for zo > 100 
and a final monotonic increase towards 1 at large r. The function gi2{r), has a peak near 
r = (Tj, obviously due to the attraction between ions of opposite charges. For r > l.bRc and 
Zq > 100, the monotonic decrease of this peak towards 1 is modified, gi2{r) has a minimum, 
smaller than 1, at r ~ 2 — 3Rc, followed by a local maximum at r = 5Rc, and for larger r 
decreases monotonically towards 1. 

The functions g'oi('") and gQ2{r) when zq < 100 have the expected behavior at large 
r characterized, respectively, by a monotonic decrease for goi{r) and increase for go2ir) 
towards 1. But, above Zq > 100, goi{r) {go2{r)) has, a local minimum (maximum) at 
r ~ 4Rc, followed at r ~ 7Rc by a a local maximum (minimum) and at large r decreases 
(increases) towards 1. 

In summary, for low values of zq, the asymptotic behavior of the correlation functions 
seems exponential, as in the case of infinite dilution of the colloids. For a large value of zq, 
it presents oscillations, although, within the limit of numerical accuracy, only one or two 
oscillations can be observed unambiguously. 

When zi = —2 and ^2 = 1 or ^1 = — 1 and Z2 — 2 , similar analysis of the asymptotic 
behavior of the correlation functions can be made with identical conclusions, since in these 
cases, the HNC equations can be solved for zq < 260 and zq < 465, respectively (cf. Table 
1). The main remark concerns the system with zi — —2 and Z2 — 1, where the onset of a 
peak in goo{r) at r ~ 2.08i?c for zq > 250, indicates that the screening of the colloid charges 
by negative ions has a range sufficiently short for two colloids to be separated by a distance 
of one or a few ion diameters (cf. Fig. 13). 

At the packing fraction rjc — 0.0001 and — 0.001 M, the HNC equations can be solved 
numerically for symmetric {zi — —1 and Z2 — 1) and asymmetric {zi — —1 and Z2 — 2, 
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FIG. 14: At r/c = 0.0001, zi = -1, Z2 = 1, p,, = 0.001 M and zq = 460, bottom : correlation 
functions gu{r) (solid line), gi2{r) (dashed line) and g22{r) (dash-dotted line), insert : close view 
for r > 2Rc, and top : goo{f) (solid line), goi{r) (dashed line) and go2{r) (dash-dotted line), insert 
: close view for r > 5i?c- 

zi — —2 and Z2 — 1) ions up to zq < 460, zq < 630 and zq < 270, respectively. The main 
difference, compared to the case of the colloid packing fraction of rjc — 0.001, is that, for 
Zq > 100, a peak or minimum appears at r > 2Rc in the ion-ion correlation functions and 
at r > 4i?c or r > 6Rc in the ion-colloid or colloid- colloid correlation functions. Beyond this 
peak or minimum, the correlation functions go monotonically towards 1 within the numerical 
accuracy of the HNC solutions. For the three considered cases of ionic charges, near 2;^"^, 
a peak appears in goo{r) for r ~ 2.08i?c, indicating that the screening of the colloid charges 
allows a close approach between colloids (cf. Fig. 14). 

In the other studied mixtures at ?7c = 0.01, and rjc = 0.1 with symmetric and asymmetric 
ionic charges, the correlation functions have asymptotic behaviors quahtatively similar to 
those found for lower colloid packing fractions. The main quantitative difference is that the 
non-monotonic decay of the correlations appears at low colloid charges Zo ~ 10. For the 
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colloid packing fraction rjc = 0.1, the correlation functions have an oscillatory behavior at 
large r for zq > 20. At this value of rjc and 2:0 > 70 — 80, several oscillations can be identified 
at large r with amplitudes larger than the numerical error. 

From the results obtained at a finite density of the colloids, it appears that the asymp- 
totic behavior of the correlations functions, as computed from the HNC OZ equations, is 
compatible with the hypothesis that it is determined by the poles of 1/D{k) as only mono- 
tonic or oscillatory behaviors have been found. It is worth remarking that, since the type 
of asymptotic decrease is the same for all the correlation functions, it should be sufficient 
to consider only one of the correlation functions for determining this behavior, for instance 
gooi'f')- However, as in the case of the infinite dilution colloid system, such a determination, 
due to the present numerical accuracy on goo{r) can be performed only in the range of r for 
which |5'oo('") ~ 1| > 10~®. A similar remark obviously applies to all the other correlations 
functions. With the hypothesis that the asymptotic behavior of the functions Qapij) is de- 
termined by simple poles near the real axis, and considering only the two nearest poles to 



the real axis, a function gai3{r) writes 



gap{r) ^ l + A^^exp(-fi;pr)cos(fcJr + d^^)/r 

+ AlpeM-^or) cos(/c°r + <^)/r . (17) 

The relative values of Hp and determine the dominant contribution to the asymptotic 
behavior. When, for instance, Kp « Kq, and and c?^^ equal to or different from zero, 
the decrease at large r is exponential or damped oscillatory ; when obviously 
the type of decrease can be difficult to identify from the numerical values of Qaisi^), in 
particular when A^^ « Figure 15 shows the behavior of the functions gai3ij) for the 
suspensions with rjc = 0.1, zq = 108, zi = —1, Z2 = I and = 0.003M. All functions present 
damped oscillations for r > QRc, in particular, goo{r) which is above the level of numerical 
uncertainties until r ~ 15i?c- This result indicates that, assuming Hp < and a value of 
Aqq estimated from Eq. (fT^ . the asymptotic behavior of goo{r) is not determined by the 
Debye screening since for r ~ 12Rc, the factor exp(— K^r) is equal to ~ 10^^^. In Fig. 15, 
as expected, the oscillations in all correlation functions have an identical period, although 
the numerical accuracy and short range contributions for r ~ 4 — 5Rc can preclude an exact 
coincidence of the oscillations which, however, is realized with an accuracy of 10"'^ -Rc for 
r > 5Rc. 
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FIG. 15: At Tjc = 0.1, zi = —1, Z2 = I, Pv = 0.003 M and zq = 108, bottom : correlation functions 
guif) (solid line), gi2{'r') (dashed line) and g22{r) (dash-dotted line), middle : gQo{r) (solid line), 
gQ\{r) (dashed line) and go2{f) (dash-dotted line), top : for r > 5Rc, close view showing the identity 
of the oscillation period of the six correlation functions, after an appropriate shift on r (cf. Fig. 
11) of the ion-ion (dashed line) and colloid- ion (solid line) correlation functions, (notice that the 
maxima and minima of the oscillations of goo{r), goi{r) and gQ2{i") can be outside the limits of the 
figure) . 



Another example where the asymptotic behavior does not seem to be determined by the 
Debye screening, is found for rjc — 0.001, Zq — 300, Zi — —1, Z2 — 1 and p„ = 0.001 M. 
This conclusion is obtained from the fact that g'oo(^) is significantly larger than 10""^ for 
r ~ 10 — 12i?c, a domain of r where exp(— K^^r) ~ lO"^'^ (cf. Table 1). For this suspension, 
the asymptotic behavior of the correlations seems well described by a superposition of two 
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exponentials {k^ = 0, = 0) where and A°^p have opposite signs (cf. Fig. 12). 

At Tjc = 0.0001, zq = 460, zi = —1, Z2 = I and = O.OOIM, the asymptotic behavior 
is described by one exponential with a screening length Kp compatible with the value of the 
Debye length, because the maximum of goo{r) exists at r ^ 8Rc where exp(— K/^r) ^ 10^^. 
However it is worth remarking that the signs of the A^^ coefficients for the colloid-colloid 
and ion-colloid functions are opposite to their expected values ; for instance Aqq is positive 
(cf. Fig. 14). Clearly, this fact seems to preclude the interpretation of A^q in terms of an 
effective charge which implies A^q oc — Z'^ff (cf. Eq. (0)). 

For the considered suspensions, the OZ HNC equations predict that the asymptotic cor- 
relation of the two-body functions has several possible behaviors. For zq > 100, the decrease 
of these correlations at large r generally does not correspond to that induced by the Debye 
screening. Due to the large difference of size between colloids and ions, the colloid-colloid 
correlations, which have a range of several colloid diameters, determine at large r the behav- 
ior of all the correlations between the ions and the colloids and ions, as a consequence of the 
exact OZ relations. In particular, when zq increases, the excluded volume effects between 
the colloids are reinforced and gooi'f') presents oscillations and looks like the pair correlation 
of a fluid of soft core particles at low or moderate densities. 



IV. CONCLUSION 



At the infinite dilution of the colloids, the ion-colloid pair correlation functions computed 
from the OZ HNC equations allows to define effective colloid-ion potentials and colloid-ion 
charges (Zg//). The results obtained for ionic solutions where the ion density varies between 
O.OOIM to O.IM are in fair agreement with the effective potentials and charges obtained from 
PB theory when zq is smaller than 1000. For larger zq, the PB and HNC values of Zgff differ. 
While the PB values go towards a constant value both for symmetric and asymmetric ionic 
solutions, the HNC values decrease and, in the case of the asymmetric solution zi = —2 and 
Z2 = 1, Z^jj'~^ becomes negative. This last result seems to be due to an overscreening of the 
colloid charge. The HNC theory shows that, when the ionic concentration is of the order 
of 0.1 M, the corrections to the Debye screening length due to finite density effects are not 
negligible in the computation of ^e//- Although the difference between and is only 
5%, it induces a factor of 3 on the value of Zefj when in Eq. p2|l is substituted to Hp. 
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At finite concentrations of tlie colloids, all the correlation functions present an identical 
asymptotic qualitative behavior. The colloid-colloid correlations determine the decrease of 
this behavior. The analysis of the asymptotic behavior does not seem to allow, as in the 
case of infinite dilution, to define in a simple way a colloid-colloid effective potential or 
colloid effective charge from f™^(r). The main reason is that, at large r, the colloid-colloid 
correlations, for large zq values may be dominated by excluded volume effects similar to 
those existing in fluids of particles interacting by a hard or soft core repulsive interaction. 
This result is in agreement with the fact that the range of the correlations is larger than that 
expected from the value of the Debye length controlling the charge screening. Obviously, it 
leaves open the possibility, largely discussed in the literature, to define an effective colloid- 



1 screening and excluded effects from a 

37, a^- 



colloid potential, able to reproduce adequately bot 

36, 



relation different from \og{goo{r)} = — /5f™o('^) 

The HNC approximation has been proved reliable to describe the coulombic fluids by 
comparison with MC simulations. For the present model, such a test of the validity of HNC 
approximation seems prevented by the difficulty of realizing an efficient MC sampling of the 
configurations of colloid suspensions where the ratio of colloid and ion sizes is 80 and the 
colloid charges are much larger than 100. However, if the HNC approximation generally 
gives good quantitative results when it can be solved, the stability limit of the numerical 
solutions can not be faithfully associated with a thermodynamic transition. Thus, at the 
colloid packing fraction rjc = 0.001, the loss of stability of the numerical HNC solutions 
corresponding to the onset of a peak in goo{r) at r ~ 2.08, cannot unambiguously be 
interpreted as the formation of colloid aggregates. 

From the present results and those published in the literature jsO, 40 1 it seems possible to 
extend the use of the HNC OZ equations for studying colloidal suspensions with colloid-ion 
size ratio larger than 80 and, also, to consider ionic solutions with a discrete solvent, at least 
for an infinite dilution of the colloids and to determine the influence of the solvent on the 
ion and solvent-colloid correlation functions and the values of Z^ff. 
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